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<^ | Abstract 

o . 

The pion decay constant f n plays a crucial role in many areas of low energy particle 
physics. Its value may e.g. be deduced from experimental data on leptonic pion de- 
cays. Here, we provide comments on several aspects of this evaluation. In particular, 
we point out that at the present level of experimental accuracy, the value of f w is 
Q_i' sensitive to the value of the pion mass chosen in its chiral expansion. 
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O 1 Introduction 

o : 

In the framework of QCD, the pion decay constant f n may be defined through 
the coupling of the axial current to the pion [1] , 

"3 ■ (0K(0)k»> = iPnf* ; A, = u W d , (1) 

where \ir~(p)) denotes a one-pion state with normalization (tc~ (p')\ir~ (p)) = 
2(27r)V ^(p'-p). 

The pion decay constant plays a crucial role in many areas of low energy 
particle physics. First of all, it dictates the strength of leptonic pion decays, 

tt~ -»■ fvi\ £ = e,fi , (2) 
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with rat^ZI 

r<°>(7T iv) = GW f fl mM (l - 4V (3) 
87T \ m*J 

in the absence of electromagnetic interactions. It also enters the Goldberger- 
Treiman relation [2], 

fn g-KN = V2m p g A , (4) 

which relates the weak and strong coupling constants /„., g& and with the 
proton mass m p . This relation is exact in the chiral limit m u = rrid = [3] 
and may well be considered the starting point of precise low energy hadron 
physics. It has been and still is a crucial test of our understanding of low 
energy hadron dynamics. Further, /„. happens to determine the strength of 
Ten interactions and thus acts as an (inverse) expansion parameter in Chiral 
Perturbation Theory (ChPT) [4,5]. Last but not least, f„ is now amenable 
to lattice calculations, see e.g. Ref. [6] for an evaluation with Nj = 2 + 1 + 1, 
and for further references. 



The most precise determination of f n is presently obtained from leptonic de- 
cays of the pion. In fact, experiments have reached a level of precision which 
makes it mandatory to include radiative corrections, and to consider the rate 
for 



vr -> £u(j) . (5) 

According to PDG [1], 

U = (130.4 ± 0.04 ± 0.2) MeV . (6) 

This value is based on data [branching fraction for ir~ — > // _ 1/^(7) and lifetime], 
and on theoretical work performed in Refs. [7,8,9,10,11,12,13,14,15]. The first 
(second) uncertainty is due to the uncertainty in the value of V uc i (to the 
uncertainties in the higher order corrections in the evaluation of the matrix 
element for the decay (5)). It is the main purpose of this Letter to comment on 
recent calculations [7,9,12,15] of the transition matrix element for this process, 
and on the value of f n reported in Eq. (6), see Section 4 for details on the 
questions investigated here. 



2 The effective Lagrangian 

A very elegant and convenient tool to perform the calculation is the effective 
field theory framework set up in Refs. [8,9,12,15]. We adhere here to this 

1 To ease notation, we often write in the following it — > £v, which stands for 7r~ — > 
l~i>£ or for ir + — > £ + f£. 
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method, and come back to its relation to the underlying theory below. The 
pertinent lowest-order effective Lagrangian for three flavours can be found in 
Ref. [12]. Here, we consider its two- flavour version, 



£e// =-4- <«X + X+) + e 2 F^Z (u^Qu 2 Qu^ 



(7) 



where the flavour-trace is indicated by (), and 

r° 

y/2n~ - k'- 



u 







(8) 



Here, F denotes the pion decay constant in the chiral limit, in a normalization 
which is standard in ChPT: Let / = fir\ mu=md= o- Then / = \/2F. In the 
absence of pseudoscalar densities, one has 



X+ = u^xd + U XU , X = 2Bdi&g(m u , m d ) . 



(9) 



u„, = 1 



u ' 



The external vector and axial vector external sources and a M contain also 
the lepton and photon fields, 

,f (dp - ir„) u-u(d lt - ilu) w f ] , 

i„ - eQA^ + (h^£LQ W + /i.e.) , 

£=e,/j. 

v^ + a^- eQAf, , (10) 



- a, 



with 



Q 




Q w = -2V2d 




(11) 



For further notation, and for the terms at next-to-leading order, see Ref. [12]. 



3 7T — > iv without electromagnetic interactions 

We first consider the non-radiative decay n — > Iv in the absence of electromag- 
netic corrections, and set e = in the effective Lagrangian. It is convenient to 
use the axial current as an interpolating field for the pion, 

out|4,(0)t|0) = ^TZ^T( qi , q 2 ) ; p = Ql + q 2 . (12) 
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At p 2 = m 2 , the quantity T(qi, q 2 ) is the transition amplitude for n — > £ up, 



As has been pointed out by the authors of Ref. [15], the chiral expansion of 
the transition amplitude, evaluated in the effective field theory framework 
mentioned, amounts to the chiral expansion of f w . A proof of the statement is 
provided in the Appendix. It relies on the fact that the pion decay constant 
also occurs in the correlator of two axial currents, 



The quantity R is holomorphic in the complex p 2 -plane, cut along the real axis 
forp 2 > 9m 2 . [The ambiguities inherent in the definition of A 11 ", generated by 
the short distance singularities in (0\TA»(x)A^(0)\0), affect R and C only] 
One could as well use the relation (15) to define f w in pure QCD. 



4 Two questions 

Because the effects of real and virtual photons have to be included for the 
confrontation of theory with experiment, infrared (IR) singularities occur in 
intermediate steps of the calculation. One source of these singularities is the 
fact that the correlator A 11 ", evaluated at e ^ 0, develops a branch point in 
the form factor B(p 2 ) at p 2 = m 2 , as a result of which there is no isolated 
pole-contribution - the decomposition Eq. (15) does not hold anymore in the 
presence of electromagnetic interactions. In addition, the formula Eq. (12) 
cannot be used without further ado when e ^ 0. So, one may wonder about 
the role of the pion decay constant that is determined in pion decays: 

i) What is its relation to the correlator A^ v at e ^ 0? 

ii) What is its relation to the pion decay constant f n in pure QCD, as it 
occurs in the decomposition Eq. (15)? 

To the best of our knowledge, these two questions were never discussed in full 
detail in the literature, and we find it instructive to shed additional light on 
the issue. On the other hand, we do not question the final algebraic result for 
the rate as provided in the works mentioned - we have nothing to add here. 

5 Photons generate a branch point 

5.1 Infrared regularization 

To perform the calculations, one may tame the IR singularities by providing 
the photon with a small mass m 7 , such that the decomposition Eq. (15) still 



T = iV^GFVudf^rriiuiq^v^) . 



(13) 



B{p 2 ) 
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g) ~ ~) 

Fig. 1. Diagrams that contribute to the correlator at one- loop order. Double 
lines denote axial currents, solid (wavy) lines pions (photons). Contributions from 
counterterms are not shown, nor do we indicate graphs which generate local terms 
only. 

holds, with f n — > f n , where the constant f n now also includes contributions 
from virtual photons. The remainder R(p 2 ) generates a branch point at p 2 = 
(m n + m 7 ) 2 . The quantity f n is gauge dependent and diverges logarithmically 
as the photon mass is sent to zero. Further, the LSZ formula Eq. (12) remains 
true, with f n — > /„.. The IR singularities cancel at the end when adding the 
rates for tt — > iv and for i\ — > lv(n^\ and sending m 7 to zero at the very 
end of the calculation provides the desired result. This is the method used in 
Refs. [9,12,15]. 

Because dimensional regularization is a very useful ultraviolet and infrared 
regulator for our purpose, we adhere in the following to this alternative reg- 
ularization [16,17,18,19], where the photon mass is set to zero from the very 
beginning. We start the discussion with the evaluation of the correlator 

^r=i/^<°™^(°)i°>« us) 

in d space-time dimensions, in the presence of electromagnetic interactions, in 
the framework of the effective theory defined by the Lagrangian Eq. (7). The 
index e indicates that virtual photons are included. 



5.2 Loop contributions at d ^ 4: explicit expressions 

To keep everything as simple as possible, we restrict ourself to a one-loop cal- 
culation. Some of the pertinent graphs which contribute to A^ u are displayed 
in Fig. 1. A typical contribution generated e.g. by diagrams Figs, la+le reads 
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where J(p 2 ) denotes the one-loop integral 
d d l 1 



Ap 2 



1 r d a l 
~iJ (2-kY 



(2n) d (M 2 - (p-l) 2 )(-l 2 ) ' 



(18) 



As for the notation, we recall that the chiral expansion of the charged pion 
mass starts out with M 2 + = M 2 + 0(e 2 ,m 2 ), where M 2 = (m u + m d )B, see 
Refs. [8,9]. The ellipsis in Eq. (17) stands for additional contributions. The 
double pole at p 2 = M 2 is removed by mass renormalization in the standard 
manner. We split off the part which contributes to the pion mass, 



J(p 2 ) = J(M 2 ) + (1 - z) J(p 2 ) ; z 
such that 



M 2 ' 



Ip^PvF 2 
Ml - p 2 



[1 + e 2 J(p 2 )) + ■ ■ ■ ; M 2 = M 2 (l - e 2 J(M 2 )) . 



(19) 



(20) 



If the photon were massive, the photon propagator in Eq. (18) would be 
replaced by l/(— I 2 ) — > l/{m 2 f — I 2 ), and the corresponding loop function 



J_(p 2 , m 7 ) would be of the form J_{p 2 



C + 0(p 2 - M 2 ), with C a finite 



constant at d — 4, and the decomposition Eq. (15) would hold for A^ u as well. 
However, massless photons render loop contributions singular at threshold. To 
investigate the structure of A^", we evaluate loop integrals at d ^ A,p 2 ^ M 2 . 
In particular, the quantity J{p 2 ) can be expressed in terms of hypergeometric 
functions [20]. Expanding the parts which are regular at p 2 = M 2 , we find 

J(p 2 ) = M 2w f{w) (l - - z) [(1 - z) 2w 9l (w, z) + g 2 (w, z) 



92 



T{-w) 



, 9i 



r(i + w ) 



1 + 2w (4vr) 2 +^ 
\- T(n-w) _ _! 



yl+W 



d 



-gi + w h(d, z) , w = '- - 2 



(21) 



The function h(d, z) is regular at d = 4, z = 1, and J_(p 2 ) behaves as 
I(P 2 ) 



L = 



(L + k) [(1 - z) d - 4 - 1 

d-4 r ! 1 



+ 



N 



d-A 



-(ln4vr + r'(l) + l) 







k- J- 




' 27V 


/i 2 



AT = 16tt 2 



(22) 



in the vicinity of p 2 = M 2 , d = 4. The renormalization scale is denoted by /i, 
and the ellipsis stands for terms that are irrelevant in the following. The region 
where this representation is valid is indicated in Fig. 2a with the shaded circle 
[use M 2 + = M 2 ]. 
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Fig. 2. The one-loop integral J and the correlator A^f in the d,p 2 plane. In Fig. a), 
the shaded circle indicates the region where the expression Eq. (22) is valid. In 
Fig. b), we indicate 3 different limiting procedures. The limit 1 (2) is relevant for 
renormalization (for the LSZ formalism). In the limit 3, J and Ae U do not exist. 
For J use M 2 , = M 2 . 



5. 3 Three limiting procedures 

We now consider three different limiting procedures [21]. First, we let d — > 4 
off the mass shell p 2 = M 2 . This limit is relevant for renormalization and is 
indicated by the path ® in Fig. 2b. Second, for the LSZ formalism, one goes 
to the mass shell first and then considers d — > 4 + [path (2) in the figure] , see 
Refs. [18,19,22]. The result is 

, [W^zyHnn -z) [limit®] 

J(p 2 ) = {" ' 1 ' [ WJ (23) 

-(L + k) [limit©]. 



In the limit (3) , J does not exist. 

As a result of Eq. (23), the correlator generates a branch point at p 2 = Mf, 



2VilVvF \ e2 +■■■ [limit©]. (24) 



2(^ _ pL^i 167r 



On the other hand, going on the mass shell p = M{ at d > 4 results in 

K = ^J(l-e 2 (^ + fc) + ---) [limit©]. (25) 

In other words, the standard renormalized contribution generates a branch 
point in the propagator. Going to the mass shell at d > 4 results in the 
standard pole behaviour of the propagator. 
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6 The correlator A% v at one loop 

Including all the graphs at one-loop order results in the following expression 
for the correlator in the first limit, 

= pV4 [limit (1)], 

7T + 

2F 2 (l + ai + e% + 0[M 4 , e\ M 2 e 2 }) , 



«i = 



1 /»,2, M2 ,,2 , M ^+\ 2M \r 



2 



6i = ^[-3 + ^(1 + In ^)] + ^ r (/i), 

1 2D 
gi = jj(6- 20 , = Y (k[ + fc£) + 4*5 . (26) 

Here, (&[) are LECs in the effective Lagrangian at order p 4 (e 2 p 2 ), see 
Refs. [5] ([23,24]). We have used the result Ref. [25, Eq. (CIO)] for the renor- 
malization of kg in any gauge. The above expression reveals the branch point 
at p 2 = M 2 + , with strengths <?i,fti which are gauge dependent. The ellipsis 
indicates terms that are less singular at p 2 = M 2 + , and terms proportional to 
g^ v . The branch point is generated by Fig. le alone - this explains the fact 
that a branch point with identical strength gi is also present in the electron 
propagator in the framework of QED [22], or in the two-point function of two 
charged pion fields in scalar QED [8]. In Ref. [22], it is shown that for the 
electron propagator, the leading infrared singularity at any order in e is given 
by the coefficient g\. 



Performing the second limit gives 

P IM P V K 2 

M 2 + - p 



k 2 



(6 - 20L + _( 3 In — - 4) + K r (^) . (27) 



= 2F 2 (l + ai + e 2 b + 0[M\ e 4 , M 2 e 2 }) , 
1 , , M 2 

The ellipsis has the same meaning as before. It is seen that in this case, the 
correlator does have - aside from less singular terms - a pole contribution at 
p 2 = M 2 + , with a residue that is divergent at d — 4, and gauge dependent. 
At e = 0, the quantity k — K\ coincides with the pion decay constant f n , 
evaluated at one-loop order in ChPT [5], 

(1 / -' \ /-' /Lf2 \ 

1 -16^ ln ^ + ^ l4 + 0{M4) ) ■ (28) 



Neufeld and Rupertsberger [9] have evaluated the pion decay constant with 
the photon mass m 7 as an infrared regulator. With the identification 
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1 TP 2 

£ --2iV< k # + 1) ' <29) 

the quantity k/v^2, evaluated at £ = 1, agrees with their F n ±, translated to 
the SU(2) x SU (2) case by use of the matching relations for the LECs l\, k r 12 g 
worked out in Refs. [26,27]. 

In Refs. [11,28], it is shown that the LEC kg depends logarithmically on the 
scale of the underlying theory [we use the matching of kg to K[ 2 as worked 
out in Ref. [27]], see also Ref. [14]. In other words, the residues k , are 
not uniquely defined in the framework of QCD+QED. An analogous scale 
dependence of kg occurs in the effective theory of the linear sigma model, 
coupled to electromagnetism [25]. 



7 Width for tt ->■ &/( 7 ) 

We come back to the evaluation of the width for leptonic pion decays (5). The 
reduction formula reads 

(P<(<Z 2 )r (gi); out|A M (0) t |0) e = M 2 KP _\ T e (q u q 2 )-p = q l + q 2 (30) 

for the non-radiative part, and similarly for the case when a photon is emitted 
in addition. Here, it is understood that the limit p 2 — > M 2 + is taken at d > 4. 
Let us denote by class I (class II) the set of graphs which do not (which 
do) contain a virtual charged lepton, and disregard the leptonic counter term 
contributions for a moment. The graphs in class I are obtained by replacing 
one of the axial currents in Fig. 1 by ^ - 7 M p££,. From the previous discussion 
and from Eq. (13), it is clear that the graphs in class I generate the amplitude 

i\/2G F V ud Km t u(qi)v L (q 2 ) = CjK , (31) 

see also the Appendix. The graphs from class II are displayed in Fig. 3. These 
contributions are of order F in this order of the momentum expansion. This 
remains true including the leptonic counterterms, and we denote the sum by 
Cue 2 F. Finally, one has to add the effect from real photon emission - the 
pertinent amplitude is denoted by CmtF . All in all, the width at one loop is 
obtained from 

r( vr _> £jy ( 7 )) = ((\ Ci k + C n e 2 F\ 2 )) + ((\C HI eF\ 2 )) , (32) 

where the symbol (( )) denotes phase space integrations, including all kinematic 
factors. The IR singularities and the gauge dependence cancel out in r(7r — > 
iv (7)), and one ends up with the expression first given in Eq. (5.1) of Ref. [12], 
matched to the SU (2) x SU (2) case considered here. Further, as pointed out 
in Ref. [14], the above mentioned scale dependence of kg is cancelled by the 
scale dependence of the leptonic LEC Xg introduced in Ref. [12]. 
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Fig. 3. Graphs which contribute to T e and which contain a virtual lepton. Double 
(solid) lines: axial current (charged pion). Long dashed (short dashed) lines: charged 
(neutral) lepton. Wavy lines: photon. 

8 The answers to the questions in Section 4 

We are now prepared to answer the questions raised in Section 4. 

Answer to question i): 

The relation of the pion decay constant to the correlator of two axial currents 
is evident from Eq. (32): the first term on the right-hand side in Eq. (32) 
contains the residue of A^ v at d > 4, which itself contains - among other 
contributions - the pion decay constant f n in pure QCD. At this order in the 
low energy expansion, one may therefore factor out f 2 in the expression for the 
width [12]. The final result may be written [12] in the form given by Marciano 
and Sirlin (Eq. (7a) of Ref. [7]), 



Gl\V ud \ 2 Jl_M^w: 
8tt 



r(vr _> t v ( 7 )) = "prudi _ Zi) 2 x 




l0g ^f-| l0g ^ L + F ^)-^i 
m z p 2 



(33) 



with zi = m 2 l /M 2 +. The function F(x) is given in Eq. (7b) of Ref. [7], and the 
parameters 6*2,3 introduced there do not occur at this order in the low energy 
expansion [12]. The constant C\ can be expressed in terms of LECs and mass 
logarithms, see Eq. (5.11) of Ref. [12] [its adaption to the SU{2) x SU{2) case 
considered here is straightforward]. It is now clear that, once G^, V u d and C\ 
are known, the pion decay constant is fixed through data on tt — > Ivi^i). This 
is how the numerical result Eq. (6) was obtained by PDG [where higher order 
corrections in the width, as worked out in Ref. [15], were taken into account 
as well.]. 

Answer to question ii): 

Here arise two subtle points. 

a) In factoring out f n in Eq. (33), one makes use of the chiral expansion of 
f n as given in Eq. (28). The choice of the mass M in that expansion amounts 
to a convention: one may use either the neutral or the charged pion mass. 
Scrutinizing the calculations performed in Refs. [12,15], we find that the value 
Eq. (6) corresponds to the choice M = M^o ~ 135 MeV [physical value of the 
neutral pion mass]. Factoring out f n evaluated at M = M n + = 139.57 MeV 
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[physical value of the charged pion mass] amounts to a renormalization of the 
constant G\. The formula (28) allows one to determine the difference between 
the two cases, 



f n (M n+ ) = f n (M 7T o) + 



h = - In 




M 2 

7T 



.0 



(U - l) + 0[p\ (M 2 + - M 2 ) 2 ] 



(34) 



Using the value Z 4 = 4.4 [29], we find 



/ ff (M w+ ) = / w (M ff o) + 0.4MeV. 



(35) 



The difference is thus quite significant - about twice the uncertainty reported 
in the PDG- value Eq. (6). The induced change in C\ is 



b) Concerning the second point, we note that, extracting a value for f n [defined 
in pure QCD, e.g. through Eq. (15)] from leptonic pion decays, where real and 
virtual photons are included, requires that one performs a splitting between 
strong and electromagnetic effects. This splitting is ambiguous [30] - the result 
depends on the procedure chosen. It is at this stage that the matching of the 
effective theory to the underlying theory matters. For a detailed analysis of this 
fact in a case which can be analyzed in a perturbative manner [linear sigma 
model, coupled to electromagnetism] , we refer the reader to Ref . [25] . Here we 
note that, in the language of the frameworks used in Refs. [11,12,13,14,15], 
the quantity depends on the scale of the underlying theory [31,25]: fJ-o^ = 
0[e 2 m q ]. A different method to perform the matching consists in evaluating 
quantities in QCD with values of the parameters relevant in QCD+QED at a 
scale fix- m this case, strong quantities are scale independent, but do depend 
on the matching scale /ii. This scenario is discussed in detail in Refs. [25,32], 
see also Refs. [33,34]. For the pion decay constant in the chiral limit, the 
dependence on fii can be worked out in the framework of QCD, using the 
observation that / is proportional to the renormalization group invariant scale 
of QCD, and applying the formula Eq. (11.6) in Ref. [30]^J. It turns out that, 
if \i\ is changed by a factor of 2, / changes by a negligible amount of about 8 
keV. The effect is so small, because the electromagnetic renormalization of the 
strong coupling constant g starts out at two-loop order and is 0(e 2 g 3 ) [30]- 
the one-loop contribution, which would be of order e 2 g, vanishes. In the linear 
sigma model, the scale dependence of / is more than one order of magnitude 
larger [25]. 

2 We are indebted to H. Leutwyler for pointing this out to us. 



d d + 2.8 . 



(36) 
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In our opinion, it would be very useful to perform the matching of the effective 
theory constructed in Refs. [8,9,12,15] to the underlying theory [the Standard 
Model] in this setting, which is used in [34]. 



9 Effects from m d ^ m u 

Finally, we comment on m d — m u effects in pure QCD, and note that the 
difference m,i — m u can occur only with even powers in f n [35]. One-loop con- 
tributions in ChPT are linear in the quark masses, and thus cannot contain 
isospin breaking terms, whereas they do occur at and beyond two-loop or- 
der. Indeed, the latter have been evaluated in Ref. [36] and found to be tiny, 
fn+ffno — 1 ~ 0.7 x 1CT 4 . Barring unexpectedly large higher order contribu- 
tions, one concludes [36] that f n + ~ f n o in pure QCD, to a very good approx- 
imation. A numerical estimate of the contributions at order e 2 in /„.+ / f n o — 1 
can be found in Ref. [37]. 



10 Summary 

We have evaluated the correlator of two charged axial currents at one loop in 
ChPT, including virtual photons, for any value of the gauge fixing parameter 
£. As is seen from the result Eq. (26), photon loops modify the holomorphic 
properties of the correlator in a fundamental manner: the pole at the charged 
pion mass becomes a branch point, with strengths «i,<?i that are gauge de- 
pendent. The residue k\ furthermore depends on the scale of the underlying 
theory [QCD+QED], through the LEC k r 9 [11,14]. For the evaluation of lep- 
tonic pion decays in the framework of the effective field theory framework set 
up in Refs. [8,9,12,15], one may use the standard LSZ formalism also at van- 
ishing photon mass, provided that the graphs are worked out in d space-time 
dimensions, and provided that the residue is evaluated at d > 4 [18,19]. We 
have discussed the manner in which the pion decay constant shows up in the 
final formula for the decay width, and have pointed out that the value reported 
in Eq. (6) corresponds to the case where the pion mass is identified with the 
neutral one. Evaluating f w at the charged pion mass increases its value by 
about 0.4 MeV. Furthermore, we note that the value Eq. (6) is based on a 
matching procedure which differs from the one used in Ref. [34]. In particular, 
f n depends on the scale of the underlying theory. 
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Appendix 

Here we prove Eq. (13), which is true in the absence of electromagnetic interac- 
tions. In this case, there are no virtual leptons in the corresponding Feynman 
diagrams, and one may evaluate the transition amplitude by considering the 
lepton current as a classical external field which we denote by X M , such that 

In = - a M - , r M = Vf, + a M , (Al) 

and 

4// = ^ + X+> (A2) 

is the pertinent leading order Lagrangian, which has the structure 

Ceff = C nn + (£^) + faR?) + 0[r\ l\ rl) ; C n7r = C eff{l=r=0 , (A3) 

with operator-valued currents L^R^. Let L M ± = 0±, with 0±(</>) = 
±0±(— 4>). At f M = 0, the terms quadratic in the external fields in the S- 
operator are 

S = J dxdyTe l J c ^ dz (a,0») x ((a p + X p )Ot) y , (A4) 

up to contact terms, which do not contribute to the matrix element in ques- 
tion. The symbol T denotes time ordering. From Eq. (A4), it is seen that the 
coefficient of axial current x axial current is identical to the coefficient of ax- 
ial current x lepton current. It is furthermore clear from the derivation that 
this statement remains true in the presence of additional terms in the La- 
grangian which contain lepton fields exclusively through their presence in the 
left current Z M . In particular, the coefficients are still the same in the presence 
of strong counterterms. From these observations and from Eqs. (12,14,15) 
follows Eq. (13). Eq. (31) is true, because the above statements also hold 
in the presence of an external electromagnetic field and of the mass term 
e 2 F A Z (u^Qu 2 Qu^). 
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